THE SPECTRAL SHIFT FUNCTION AND LEVINSON'S THEOREM 
FOR QUANTUM STAR GRAPHS 
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Abstract. We consider the Schrodinger operator on a star shaped graph with n edges 
joined at a single vertex. We derive an expression for the trace of the difference of the 
perturbed and unperturbed resolvent in terms of a Wronskian. This leads to representa- 
tions for the perturbation determinant and the spectral shift function, and to an analog of 
Levinson's formula. 



1. Introduction and main results 

1.1. Introduction. This article focuses on the study of the spectral shift function and a 
Levinson theorem for Schrodinger operators on star shaped graphs. Quantum mechanics 
on graphs has a long history in physics and physical chemistry |21,22j, but recent progress 
in experimental solid state physics has renewed attention on them as idealized models for 
thin domains. A large literature on the subject has arisen and we refer, for instance, to the 
bibliography given in [3,11]. 

A star graph is a metric graph T with a single vertex in which a finite number n > 2 of 
edges ej are joined. We assume throughout that all edges tj are infinite and we identify 
ej = [0,oo). We assume that the potential V is a real- valued function on T satisfying 



/ |Vj(xj)| dxj < oo for all 1 < j < n, 

J e, 



where we denoted the restriction of V to the edge ej by Vj(xj) = V(x)\ t 
condition, we can define the Schrodinger operator 

Hip := -ip" + Vip 

with continuity and Kirchhoff vertex conditions 



(1.1) 
Under this 

(1.2) 



V>i(0) = ...=Vn(0)=:^(0), J») = 0, (1.3) 

i=i 

as a self-adjoint operator in the Hilbert space Li(T) = ®™ =1 ^2(ej). In (|1 .3j) we denoted 
by ipj the restriction of ip to the edge ej . More precisely, we define the operator H via the 
closed quadratic form 

h[<f>] ■= J \cjj \x)\ 2 dx + J V{x)\(t>{x)\ 2 dx, 

with form domain d(h) = H 1 ^) consisting of all continuous functions on T such that 
(pi S H (ej) for every j. If V is sufficiently regular in a neighborhood of the vertex, then 
functions (j) 111 the- operator domain of H satisfy the Kirchhoff vertex condition in (]l,3p : 
otherwise this condition has to be interpreted in a generalized sense. 

l 
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Our two main results are formulas for the spectral shift function and the perturbation 
determinant of H with respect to the unperturbed operator Hq (which is defined similarly 
as H, but with V = 0) and an analog of Levinson's theorem. Special attention will be paid 
to the existence or absence of zero energy resonances. 

There are several motivations for this study. The first one is the scattering theory of 
quantum graphs. While star graphs are certainly very special graphs, it is generally believed 
that they are a correct model example for a scattering process in the presence of a vertex. 
The direct and indirect scattering theory on star graphs has been studied in great detail 
in |12| and |15| . Our results complement theirs and, in contrast to them, we advertise 
a more operator theoretic approach including, for instance, Fredholm determinants, trace 
class estimates and Krem's resolvent formula. 

A second motivation is a line of thought that goes back at least to Jost and Pais |17] ; see 
also |23| . |13| and |20j . In these works, a perturbation determinant, which is a Fredholm 
determinant in an infinite dimensional space, is shown to be equal to a much simpler de- 
terminant, typically in a finite dimensional space, such as a Wronski determinant. While 
such formulas appear in different set-ups, there seems to be no general method of knowing 
in advance the form of the 'simpler determinant'. One of the achievements of this paper is 
to derive a new formula of this kind for a star graph. 

A third motivation comes from the general interest in zero energy resonances because of 
their key role in several diverse problems of mathematical physics; for instance, the Efimov 
effect in many-body quantum mechanics [9], the time decay of wave functions |16j and 
the convergence of 'thick quantum graphs' [13], to name just a few. We also refer to the 
survey (3|5]. In particular, we hope that our results will allow us to remove the non-resonance 
assumption in the recent dispersive estimates on star graphs [19]; see also [23] for similar 
bounds in the whole line case. 

Finally, we note that the derivation of a Levinson theorem for a graph with a finite number 
of unbounded edges was mentioned as an open problem in [7] (who considered the discrete 
case). While the compact part of the graph still has to be better understood, our analysis 
explains how to deal with several unbounded edges and will be useful, we believe, in further 
developments in this direction. 

1.2. Main results. To state our main result, namely a trace formula for the operator 
(jl.2p with vertex condition (jl.3p . we need some notations. By Hn j we denote the half- 
line Schrodinger operator with potential Vj = V\ ej and Dirichlet boundary condition at the 
origin. The self-adjoint operator 



where the form domain is given by d(fiD,j) = H Q,l {ej) = {(fij G H 1 (ej) : 4>j(0) = 0}. If the 
condition (jl.ip is satisfied, then the equation 

-u" + Vu = zu, z = C 2 



H D,j 




on L2(ej) is associated with the quadratic form 
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has two particular solutions, the regular solution ipj and the Jost solution 9j. The first one 
is characterized by the conditions 

^•(0,0 = 0, 4-(o,c) = i 

and the latter one by the asymptotics 9j(x, () = e lx ^(l + o(l)) as |C| — > oo. Both solutions 
are unique, see for instance [26]. The Jost function iOj(C) is defined as the Wronskian of the 
regular solution and the Jost solution and turns out to be Wj(Q = 9j(0, ()■ 




Figure 1. star graph T 

Our first main result is 

Theorem 1.1. Let T be a star shaped graph and assume that (II. ip is satisfied for 1 < j < n. 
Then, for the Schrddinger operator (jl.2p on Li (r) with Kirchhoff vertex condition (ll.3p , the 
following trace formula holds, 

Tr ((H - C 2 )- 1 ~{H- C 2 )- 1 ) = ~ In ( ^ f[ wtf) j , Im C > 0, (1.4) 

where K(() = £? =1 ^(0, 0/^(0, C) and Wj (() = 9,(0, Q. 

Remark 1.2. We note that identity (|1.4|) is equivalent to the identity 



Tr( ( F -c 2 r 1 -(^-c 2 r 1 ) = ^E 



" 1^(0 + %K(0 I 



™i(0 K(C) C 
which should be compared with the classical result [3[T7], see also [231126]. 

Tr ((H Dtj , - C 2 )- 1 ~ (Hn tj - C 2 )' 1 ) = |^y- (1.5) 

From equation ()1.4|) . we conclude in Section [3] an explicit expression for the perturbation 
determinant D(z) and the spectral shift function ^(X;H,Hq). We recall that the spectral 
shift function can be characterized by the formula 

/oo 
a\;H,H )f'(\)d\, 
-oo 

for any / G Cq°(R) (together with the condition £(\;H,H ) = for A < inf o-(H)). An ex- 
tension of this formula for a broader class of functions, as well as several equivalent definitions 
are discussed in Section 3. 
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In Section S] we study the low-energy asymptotics of D(z) as \z\ —> 0. This allows us to 
prove an analog of Levinson's formula for the star graph. We say that the operator H on 
L2(r), given in (|1.2j) . has a resonance at £ = if the equation —u" + Vu = has a non- 
trivial bounded solution satisfying the continuity and Kirchhoff conditions. By definition, 
the multiplicity of the resonance is the dimension of the corresponding solution space. 

Theorem 1.3. Assume that 

I {l + x)\Vj(x)\dx <oo for all I < j < n, (1.6) 

J ej 

is satisfied and, if £ = is a resonance of multiplicity one, assume that 

(l + x 2 )\Vj(x)\dx < oo for all 1 < j < n. (1.7) 



Then, 

lim ^(A) = -fiV + ^fl), (l.N) 



A^0+ \ 2 

where N is the number of negative eigenvalues of H and where m > 1 is the multiplicity if 
£ = is a resonance and m = if £ = is no£ a resonance. 

Remark 1.4. We know from Bargmann's bound that iV < oo if (II. 6ft is satisfied, [2J. We 
also know that lim^o- £(A) = — N, which is an easy consequence of the definition of the 
spectral shift function. 

2. A Trace formula for Star Graphs 

In this section, our goal is to prove a trace formula for star graphs. More precisely, we will 
find an expression for Tr (R(z) — Rq(z)) in terms of the Jost solutions Qj on the edges ej. Here 
and in the following we write R(z) = {H—z)~ l and Ro(z) = (Hq — z) -1 for the perturbed and 
unperturbed resolvent, respectively. When deriving an expression for the resolvent R(z), we 
will make use of Krem's formula for which we refer to [I] and, in particular, to an article by 
Exner |10j where this formula was used in a similar context. Thereby, we need to decouple 
the operator H which we achieve by imposing Dirichlet vertex conditions on each edge ej, 
i.e., ipjifS) = for all 1 < j < n. Then the operator (11 .2[) is decoupled and the half-lines are 
disconnected. We denote the decoupled operator by 

n 

Hoc = Hdj 

3=1 

and its resolvent by Roo(z) = (H^ — z)~ l . In what follows, we will skip for simplicity the 
indices at the coordinates and use the notation ipj{x) := ipj(xj), 1 < j ; < n, for a function 
defined on the edge ej of T. 

Proof of Theorem \l.l\ It is a well-known fact, see e.g. |26j . that under assumption (II. ip for 
all z = ( 2 such that Imz / and Wj(£) / 0, the resolvent Rdj(z) = (Hdj — z)' 1 is an 
integral operator with kernel 

R DA x,y; 2 ) ;= ^0ff(y,C) x < yi c = zl ft 
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and RD t j(x,y; z) = RDj(y,x;z). Hence, the resolvent Roc(z) is a matrix integral operator 
with the kernel 

R j?e( x >y, z ) : = 5 jt eRD,j(x,y,z), l<j,£<n. 

Having the resolvent kernel Rf p of the decoupled operator H^, we can use Krem's formula 
[l] to determine the kernel of the resolvent R(z). Let p{H) be the resolvent set of the operator 
H and p(Hq) the resolvent set of Hq. The formula states that for any £, such that Im£ > 
and z = £ 2 G p(i?oo) n p(H), the resolvent R(z) is a matrix integral operator with kernel 
Rj t i(x,y; z) = i?°^(x,y;z) + \jt0j(x,Q6e(y,Q- Li order to determine the coefficients \jg 
we proceed as follows. For any / = (fi(x), . . . , f n (x)) T £ ^(r), the function ip(x) := 
f R(x,y; z)f{y) dy has to satisfy the equation Hip = C^ip + / and the Kirchhoff vertex 
condition. This leads to a system of n linear equations for the coefficients \j£. It turns 
out that \ji = {-KlQOjfrQOtfrQ)- 1 , with K(0 = E"=i ^(0, 0/^(0, C), see also [TO]. 
Thus, 



R u {x^z) := R%{x,y;z) - (2.1) 



(x,Q0i(yX) 

K{Qe j (o,c)^(o,c) 

This representation allows us to compute Tr(R(z) — Rq{z)). First, we note that the op- 
erator R{z) — Rq{z) is a trace class operator. This can be seen as follows. As the quo- 
tient in (|2.ip is a perturbation of finite rank, we only have to show that the difference 
Roo( z ) — Roo(z) is trace class. Here R^(z) is the resolvent of the unperturbed decoupled 
operator H&P = ©j=i (—d 2 /dx 2 ) on ©™ =1 ^(cj). Similarly, we denote by R^-{z) the 
resolvent of the unperturbed operator H^- = —d 2 /dx 2 on L2(ej). Under condition (jl.lj) 
the operator ^y\Vj\ ^R^j(z)J is Hilbert-Schmidt for all a > 1/4 and all 1 < j < n, as 
can be easily checked (see e.g. [26] . Lemma 4.5.1). Hence, the Birman-Schwinger operator 
\/\Vj\ R^Qj (z) y/Vj, with \fVj := sgn(Vj)\/\Vj\, is trace class and has for z G p(Hoj) no 
eigenvalue —1. Thus, the following resolvent identity for the half-line Schrodinger operator 
holds, 



It follows from this resolvent identity that R OQ (z) — R^(z) is a trace class operator. In view 
of (|2.ip it follows that also R(z) — Rq{z) is a trace class operator and 

n „ 

Ti(R(z) - R (z)) = J2 (RD, j {x,x,z)-R { § j {x,x,z)) dx (2.2) 



j= i " «y 

' ' + — — I dx. 




j 



2 {x,Q e 2i < 



The computation of the first integral on the right-hand side is the classical Jost-Pais result 
|17| recalled in Remark 11.21 

R Dd (x,x,z) -R^ Ax,x,z)) dx = - ^|r . (2.3) 
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Here the derivative with respect to £ is denoted by a dot, "• = d/dQ n . To compute the 
second integral, we use the following equation which is true for any two arbitrary solutions 
of the equation Hojipj = C 2 "^'; namely 

2(uj(x, ()vj(x, C) = (u'j(x, C)vj(x, C) - uj(x, ()v'j(x, ())'. 

Applying this identity to Uj = Vj = Oj, we get 



9](x,C) 

ax 



9' ] (x,C)0j(x,O-O 1 (xX)O' 3 (x,O 



+ K(C)e](o,() 2CK(C)^ 2 (o,C) 

First, we consider the case of compactly supported potential Vj. Then, for large x the Jost 
solution for the half-line Schrodinger operator Hd,j is given by 9j(x, () = e^ x and we have 

efa, C) = iCe i(x , 9j(x, C) = ixe**, fyx, Q = (t - xQ)e^ x . 

Therefore, for large x, 

9 / J (x,C)9,(x,O-0 J (x,O0 , J (x,C) = ie 2 < x . 
Note that the function e 2l< * x vanishes for x — > oo, as Im£ > 0. We therefore get 

f^k + o](o,()K(0 y fr[ 2CK(c)9](o,c) 2(my { -" ] 

By density arguments (see [Prop. 4.5.3, based on the fact that ^ l /\V\{H£,fi + l)~ 1 ^ 2 is a 

Hilbert-Schmidt operator under condition (jl.ip . the result can be extended to all potentials 
Vj satisfying (jl.lj) . Similarly, we consider the case V = and obtain 

Combining (j!T3j) . and (|23j) with (pl2j) . we finally arrive at 



This is the claimed formula. 



□ 
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3. The PERTURBATION DETERMINANT and the spectral shift function 

Identity (jl.4p implies an explicit expression for the perturbation determinant 

D(z) := det(l + Wr (z)^/\V\), z G p(H ), 

where yV = (sgn V) y/\V\. Strictly speaking this is the modified perturbation determinant, 
nevertheless we shall refer to it simply as the perturbation determinant in what follows. 
Note that under the assumption (jl.ip the perturbation determinant D(z) is well-defined 
since the operator y^\V\(Ho — Z)"V2 is Hilbert-Schmidt and therefore yV Ro(z)y/\V\ is 

/ , 0) \ 1/2 

trace class. This follows as above from the fact that y \Vj\ [R D -(z)\ is Hilbert-Schmidt 
for all 1 < j < n together with f)2. 1 j) for V = as the corresponding second term on the 
right-hand side of (12. ip is of finite rank. 

Furthermore, a simple computation shows, see also (0.9.36) [25], that the perturbation 
determinant is related to the trace of the resolvent difference by 

D~ 1 (z)D'(z) = Tt(Rq(z) — R(z)), zG p(H )n P (H). (3.1) 

Hence, in view of Theorem 1 1.1 1 we conclude that 

^{z^K{z^)\[ jWj {z' 



D-\z)D'{z) 



z-i/ZKizV^UjWjiz 1 / 2 ) 



here we choose the square root of z such that Imz 1 / 2 > 0. From which it follows that 
D(z) = Cz- 1 l 2 K{z l / 2 )Y[ n j=l w j {z 1 ' 2 ), for some C G C. The coefficient C is fixed by the 
asymptotics of a perturbation determinant, namely 

lim D{z) = \. (3.2) 

I lm z\— >oo 

This asymptotics is true if the operator IVp/^-Ho — z)~ 1 / 2 is Hilbert-Schmidt, see e.g. 
(0.9.37, [25]). As |C| -»• oo, we have [Hl[26] 

^■(0 = ^(0,0 = 1 + Odd- 1 ) and K(() = ni( + 0(1). (3.3) 

This implies that C = 1/in. Thus, we have proved 

Corollary 3.1. Assume that (jl.ip is satisfied. Then, for z G p(H), the perturbation deter- 
minant of H with respect to Ho is given by 

D M=f^rW i/2 >. pw) 

where Imz 1 / 2 > 0. 

Our next goal is to determine an explicit expression for the spectral shift function £(A; H, Hq) 
for the pair of operators H,Hq in L2(T). If (jl.ip is satisfied, then by the argumentation 
above the resolvent difference R{z) — Rq{z) is trace class for all z G p{H). In this case it 
is known from general theory that for all — c < inf o~(H) there exists a real-valued function 
£ C (A) for the pair of operators R(—c),Rq(—c) such that the relation 

/oo 
£ c (A)/'(A)dA (3.5) 
-oo 
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is true for all functions / G Cq°(R). This formula goes back to Lifshits [18j . 
The spectral shift function for the pair H, Hq is defined by the relation 

£(A; H, H ) := -£ C ((A + c)" 1 , fl(-c), R (-c)) (3.6) 

for A > — c and £(A; H, Ho) := for A < — c. It can be shown that this definition is 
independent of the choice of c. By a change of variables the formula (13. 5f) for the pair R, Rq 
can then be transformed into a formula for the pair H, Hq and yields 

/oo 
C(X;H,H )f'(X)dX, (3.7) 
-oo 

for all functions / G Cg°(R). 

The next theorem among other things extends the class of admissible functions / in this 
trace formula. 

Theorem 3.2. Let H be the Schrodinger operator in L2(T) given in (11. 2ft with the Kirchhoff 
vertex condition and Hq = —d 2 /dx 2 the corresponding unperturbed operator. Assume that 
condition (jl.ip is satisfied. Then, the spectral shift function for the pair of operators H, Hq 
is given by 

£(\;H,H ) = 7T- 1 lim arg D(\ + ie), 

where arg D{z) = Imln D(z) is defined via In D(z) — > as dist(z, ct(Hq)) — > oo. 
Moreover, 

/oo 
e(A;H,H )(A-z)- 1 dA, z G p(Hq) n p(H), (3.8) 
-oo 

and fl3.Tf) holds provided f has two locally bounded derivatives and for any e > 0, m > — 1/2 
as A — > oo, 

/'(A) = 0(A _m_1 " E ), /"(A" 2 " 1 - 2 ). (3.9) 

Finally for m > —1/2, 

/oo 
|£(A;HH )|(l + |A|r m - 1 dA<oo. 
-oo 

In the proof of Theorem 13.21 we make use of some results from abstract scattering theory 
which we collect for the reader's convenience in the following proposition. These results can 
be found e.g. in [Chapter 0.9, [26] ]. 

Proposition 3.3. Let h and ho be lower semi-bounded operators on a Hilbert space and 
v = h — ho. Assume that the operator \/\v\ (ro(— c)) 1 ^ 2 is Hilbert- Schmidt for some c < 
inf (cr (h) U c(/io)) > where ro(— c) = (ho + c)^ 1 and r(—c) = (h + c)" 1 . Moreover, assume that 

/oo 
\\r(-t) - ro(-t)\\ir m dt < oo (3.10) 

for some m G (—1,0). Let the spectral shift function £(\;h,ho) be defined as above and the 
(modified) perturbation determinant by 

d(z) = det(l + sgnvy/\v\ro(z)^/\v\), z G p(h ). 

Then, d(z) — > 1 as dist(z, o~(ho)) —> oo and 

£(A; h, ho) = ir~^ lim &rgd(\ + ie), (3.11) 
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where argd(z) = Im\nd(z) is defined via lnd(z) — > as dist(z, o~(ho)) — > oo. Moreover, 

roo 

\nd{z) = £(\;h,h )(\ - z)' 1 dX, z £ p(h ) n p{h) , (3.12) 



and 

Tr(f(h)-f(h ))= I a^h,h )f'(X)dX (3.13) 



provided f is as in (|3.9p . Finally, 

|£(A;/i, /io)|(l + |A|)- m - 1 dA<oo. 



In what follows, we prove that for the case of star graphs T the condition (|3. 10[) is satisfied 
for all m > —1/2. Then, by Proposition 13.31 the assertions of Theorem 13.21 will follow. We 
need to compute the trace norm of a rank two operator. 

Lemma 3.4. Let% be a Hilbert space andf,g G T~L. Further, assume thatlZ = (•, f)f—(-,g)g 
is an operator of rank two on %. Then, the trace norm ofTZ is given by 

ll^lli = ((ll/H 2 + N| 2 ) 2 -4|(/, 5 )| 2 ) 1/2 . (3.14) 

Proof. We may assume that / 7^ 0, for otherwise formula (|3.14j) is obvious. We construct an 
orthonormal bases on Ran(7£) by applying the Gram Schmidt process, 

B = {//ll/ll, (9 - (f,g)f/\\f\\ 2 ) I (IMI 2 - llf,<?)l 2 /ll/llT 1/2 }- 

The operator 1Z is described by a 2 x 2 matrix M = (m)^, where 

mil = ll/ll 2 " l(/,9)| 2 /ll/H 2 , m 22 = -\\gf + |(/, 5 )| 2 /ll/l| 2 , 

mi 2 = mil = - ((/, <7)/||/||) (Nl 2 " l(/,9)| 2 /ll/H 2 ) 1/2 • 
The singular values of M turn out to be 

*i = ^(ii/ii 2 -ni 2 +^, *2 = ^-imi 2 + Ni 2 ), 

1 /2 

where oj := ({\\f\\ 2 + ||g|| 2 ) 2 - 4|(/,5)| 2 ) 1 . Hence, \\TZ\\i = si + s 2 = co. This proves 
(13341 . □ 

Lemma 3.5. Let H be the Schrddinger operator in £ 2 (r) given in (|1.2|) with the Kirchhoff 
vertex condition and Hq = —d 2 /dx 2 the corresponding unperturbed operator. Assume that 
condition (II. ip is satisfied. Then, the resolvents R(z) and Rq(z) satisfy for all e > and 
for t big enough, 

Proof. By adding zero we estimate the trace norm by 

\\R(-t) - i?oH)||i < H^lli + \\Roo(-t) - R&H-t)\\x, (3.15) 

where 1Z := R(—t) — i? OC) (— t) + R^(—t) — R (—t) is an operator of rank two. The second 
norm in the right-hand side of (|3.15|) can be estimated by \\Roo(—t) — Roo(— t)\\i < ct~ 3 / 2+e , 
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see [26], Lemma 4.5.6. Thus it remains to bound the norm of the rank two operator 1Z. 
From Krem's formula (|2.ip we have 

(Ro(x k ,x e ,-t) - R^(xk,X£,-t))k£ = fk(xk)fe(xe), fk(x k ) = ^^ 1/4 , 
and hence, 

Ro(-t) - R&X-t) = (-, /)/, /=(/!,..., fnf- 

Further, 

Q k (x k ,iVi) 



(R(x k ,Xi, -t) - R 00 {x k ,xi, -t)) H = 9k{x k )gi(xi), g k (x k ) 



y/-K(iy/i)O k (p,iy/i) 

and 

R(-t) - Roo(-t) = (-,5)3, g = (51, •• • ,3n) T - 
In view of (|3,14p we have 

- Roo(-t) + JS£)(-t) - Ro(-t)h = ((ll/ll 2 + Ibll 2 ) 2 " 4|(/,9)| 2 ) V2 • (3.16) 

In the remaining part we show that (||/|| 2 + \\g\\ 2 ) 2 -4|(/, c/)| 2 = 0(^ 3 ) as i -> 00. Let us set 
g = f + h and note that h is a real-valued function. Then by applying the Cauchy-Schwarz 
inequality and the arithmetic inequality, 

|2 1 ||„||2\2 A\/r „m2 Ai\\t\\2 1 It i.\MIUl2 , II z. ||4 ^ a \\ u Il2 II t\\2 , o 1 1 1, 1 1 4 



+ - mg)r = w\r + if, h ))\n z + ihi 4 < mvw + mr- (3.17) 

We compute 



E / ^—75-^ = (at)" 1 - (3.18) 
Next, we consider h = h\+h2 , where ^ = i, • • • , h k ^ n ) T ', fc = 1,2, 

^ 2 ,j = n^/-^(iVt)^(o,*vt)^ - (^ i / 4 ) _1 j e -^. 

Because of (13.31) we have as i — >• oo, 

-l 

= O^ 1 / 4 ) (3.19) 



-K(iVt)ej(0,iVt) 
and further, see e.g. Lemma 4.1.4. 



9j{x h iVt) - e'^i < -^=e-^, cel. (3.20) 

vt 

Hence, in view of (|3. 19|) and (j3.20|) . ||/ii|| = 0(i _1 ). To compute the asymptotics for ^2 we 
rewrite 

foj = (V^t 1/4 - y/-K{iy/i)6j(0,iy/ifj (J^K^)6 j {^iVt)\ f^Xj) 



n 



Vi + K(ix/i)9 2 (0,iVi 



fj (xj 



y/nt 1 ^ + J-K(iVt)9j(0,iVt) J-K{iVt)0j(0,iVi) 
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Together with (|3TTSj) and (1XT91) . this leads to \\h 2 \\ = O^" 1 ) and therefore, \\h\\ = 0(t _1 ). 
This yields in view of (|3.14p and (|3.17p that ||7£||i = 0(i~ 3 / 2 ). This proves the assertion of 
the lemma. □ 

With Lemma 13.51 all assumptions of Proposition 13. 31 are fullfilled and therefore the spectral 
shift function for the pair H, Hq in L2(T) satisfies the relations (|3.12p and A3. 13H . Especially, 
we have proved Theorem 13.21 

Remark 3.6. Lemma [331 implies that (H + c)P — (Hq + c) 13 is trace class for (3 < 1/2, see [26] . 



4. LOW-ENERGY ASYMPTOTICS AND LEVINSON'S FORMULA 

In this section we study the low-energy asymptotics of D(z) as \z\ — > 0. This will allow us 
to prove an analog of Levinson's formula for star shaped quantum graphs. Throughout this 
section we assume that (ll.6p is satisfied. 

Definition 4.1. We say that the operator H in L2(r), given in (II. 2p . has a resonance at 
C = if the equation 

_ u " + Vu = (4.1) 

has a non-trivial bounded solution satisfying the continuity and Kirchhoff conditions. By 
definition, the multiplicity of the resonance is the dimension of the corresponding solution 
space. 

Remark 4.2. We shall show below that £ = is never an eigenvalue. 

We recall some auxiliary results on half-line Schrodinger operators. For the half-line 
Schrodinger operator the Jost solution of the equation Hju = C, 2 u was characterized by its 
asymptotics 6j(x,C) = e lx ^(l + o(l)) as \(\ — > oo. Further, the function 6j(x,0) satisfies the 
equation Hju = and its behavior at ( = is given by 

6j(x,0) = 1 + (J y\Vj(y)\dy\ = l + o(l), as x oo, (4.2) 

(see e.g. Lemma 4.3.1., [26J). Recall also that the Jost function is uij{C) = 8j(0,()- If 
Wj(0) = 0, the low-energy asymptotics 

Wj (C) = -iw^C + o(C), C^O (4.3) 
(i) 

is true with some constant Wq ^ 0. 

Lemma 4.3. If Q(x,0) is a bounded solution of equation (14. lj) . then for some Cj G C ; 

n 

©M) = Ci %(x,O). (4.4) 
Moreover, the operator H cannot have a zero eigenvalue. 

Proof. If ©(x,0) solves the equation (14.1j) with £ = 0, then the restriction of @(x,0) to the 
edge ej is a solution of the corresponding zero-energy equation on the half-line for every 
1 < j < n. As stated above, the function 6j(x,0) solves the zero-energy equation and is 
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bounded at infinity by (14.20 . We note that 9j(x,0) is the only solution of the zero-energy 
equation on the half-line, which is bounded at infinity. Indeed, the solution 



for all 1 < j < n, (Lemma 4.3.2, [26]). Here xq is an arbitrary point such that 9{x) ^ 
for x > xq. Finally, equation (|4.ip cannot have a nontrivial solution belonging to -^(r) at 
infinity as @j(x, 0) = Cj + o(l) for x — > oo. 



Lemma 4.4. Let M := #{j : Wj(0) = 0}. 

(1) If £ = is not a resonance, then either M = and K(0) ^ or M = 1 and K{Q 
has a pole at £ = 0. 

(2) Assume that Q = is a resonance of multiplicity m > 1, i/ien either 

a) any resonance function vanishes at the vertex and m = M — 1 > 1, or 
6 J i/ie resonance is of multiplicity one, the corresponding resonance function is non- 
zero at the vertex and m = 1, M = and K(0) = 0. 

Proof. We first note that if M > 1, i.e. #j(0, 0) = for some j, then any resonance function 

must vanish at £ = because of the continuity condition. 

(1) If M > 2, then it is always possible to construct a zero-energy function by setting 
Cj = if Wj(0) ^ and determining the Cj's such that the Kirchhoff vertex condition 
is fullfilled if Wj(0) = 0. Hence, if £ = is not a resonance, then necessarily M < 1. 
If M = 1, then obviously K(() has a pole at £ = 0. Moreover, if M = and £ = 
is not a resonance, then i^(0) ^ 0, because if -fC(C) would vanish in £ = 0, then it 
would follow from 



that the function Q(x, 0) given in (|4.4p is a zero-energy resonance function for suitable 
Cj. Hence, K(0) ^ if C = is not a resonance. 
(2) a) If C = is a resonance, then 1, as it is not possible to construct a resonance 
function satisfying the vertex conditions and having support on only one edge of T. 
If M > 2, then because of the continuity condition any resonance function has to 
vanish at the vertex. Further, we set Cj = for all j with Wj(0) ^ 0, then there are 
M — 1 linearly independent choices for the remaining Cj's such that the Kirchhoff 
vertex condition is fullfilled. Hence, the multiplicity of the resonance function is 
m = M - 1 > 1. 

b) If £ is a resonance with M = 0, then 6j(0, 0) 7^ for all j and the coefncents Cj 
are determined uniquely by the n — 1 continuity conditions and the Kirchhoff vertex 
condition. Further, this implies because of (|4.5p that K(0) = 0. 



Jx 

which is linearly independent of 0j(x, 0) has as x — > 00, the asymptotics 

Tj(x, 0) = x + o(x) 




□ 




(4.5) 



□ 



In the following proposition, we give the low-energy asymptotics for D(z) as \z\ — > 0. 
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Proposition 4.5. Let m be the multiplicity of the resonance £ = 0, with the convention that 
m = if Q = is not a resonance. If m = 1, we assume in addition that condition (|1 .7|) is 
satisfied for all 1 < j < n. Then, as C, — ^ 0, 

D(z) = cC m -\l + o(l)), z = ( 2 , (4.6) 

with c 7^ 0. 

For the proof of Proposition 14.51 we shall need the following 

Lemma 4.6. Let Hpj = —d 2 /dx 2 + V(x) in L2([0,oo)) be given with Dirichlet boundary 
condition at the origin and assume that J*q°°(1 +x)\V(x)\dx < oo. Let 0(x,£) be the Jost 
solution on the half-line. 

(1) If 9(0,0) = 0, then 

9(0,0)6' (0,0) = -i. (4.7) 

(2) 7/0(0,0)^0 and / °°(1 + x 2 )\V(x)\ dx < oo, then 

9(0, 0)9' (0, 0) - 9'(0, 0)9(0, 0) = -i. (4.8) 

Proof. If 0(0,0) = 0, then 0(0,0) is defined for V having a first moment and 0(0,0) = 
— ico, Co 7^ 0, by [26], Proposition 4.3.7. Further, it was shown in [26J (4.3.11), that 

ip(x,0) = c 9(x,0), (4.9) 

where 93 is the regular solution of — u" + Vu = C, 2 u. Taking the derivative with respect to 
x on both sides of (I4.9P and setting x = yields that 1 = co0'(0, 0), as the regular solution 
of the Dirichlet problem was defined by the condition (p'(0,Q) = 1. Hence, equation (|4.7|) 
follows. If 0(0,0) 7^ 0, then 0(0,0) is only denned for V having a second moment and 

(p(x, 0) = i9(0, 0)9(x, 0) - i9(0, 0)9(x, 0) (4.10) 

by [26], Corollary 4.3.11. Again, taking on both sides of (I4.10P the derivative with respect 
to x and setting x = leads to ()4.8|) . 

□ 



Proof of Proposition \4-5\ We consider the explicit expression for the perturbation determi- 
nant D(z) which was given in (I3.4D . 

^)=fPn^(c), z=c 2 . (4.H) 

Let us first consider the case M = 0. Then by Lemma 14.41 either Q = is not a resonance 
and K(0) 7^ or C, = is a resonance and K(0) = 0. If Q = is not a resonance, then 
we see from (I4.1ip that as Q — > 0, D(z) = c£ -1 (l + o(l)), c 7^ 0. If ( = is a resonance, 
then we consider K(0) = £™ =1 (^(0,0)0^(0,0) - 0^(0,0)0^(0,0)) 9j 2 (0,0) which by gSJ is 

the same as K(0) = E™=i^7 2 (°>°) + °- Hence, by applying I'Hospital we have as £ — )■ 0, 
D(z)^c^0. 

Next, we consider case M > 1. Without loss of generality let 0i(0, 0) = . . . = 0j\/(0, 0) = 0. 
We rewrite (I4.11|) as 




j=l ^ ' s; fc=i j=M+l n ' s; fc=l 
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Obviously, the second term on the right-hand side is 0{Q M ). In the first term on the right- 
hand side we have for each 1 < j < M, as Q — t-0 

Ql /n A\ n n 



fc=l 
Mi 



^-^(0,0)11^(0.0) n ^(0,0) + O(C 



fe=l fc=M+l 

= C^A^rU^M II ^(0,0) + 0(C M ). 

W U )fc=l fc=M+l 

In view of (14, 7j) and (14.12P we arrive at 

y-Af-1 / A/ ™ M 

D{z) = ^-- J] 4(0,0) J] fe (O,O)f]T(0;-(O,O)) 2 + O(C) 
= cC M ~ 2 + 0{C M - 1 ), c^0, z = ( 2 . 

□ 

In the remaining part we prove an analog of Levinson's formula for star shaped graphs. 
For k G M, we set D(k 2 ) = a(k)e iri ^ k \ where a(fc) = |£>(/c 2 )|. Then it follows from the 
representation D(k 2 ) = {K{k) / {ink))Y\^ =l Wj{k) that 

-rj(k) = ?/(-*;). (4.13) 



Indeed, it follows from the uniqueness of the Jost solutions 9j(x, £) that 6*j(x, A;) = #j(x, —A;), 
and 6^(x, k) = 9'Ax, —k), and hence also Wj(k) = 9j(0, k) = Wj(—k). 

Proof of Theorem \1.3[ First, we note that there is a spectral theoretical result relating the 
zeros of the perturbation determinant to the eigenfunctions of H as follows. The function 
-D(C 2 ) has a zero in £ of order r if and only if £ 2 is an eigenvalue of multiplicity r of the 
operator H, [25]. Obviously, the zeros of D(C 2 ) lie on the positive imaginary axis as H is a 
self-adjoint operator and therefore it may have only real eigenvalues. 

We apply the argument principle to the function D{Q 2 ) and the contour T^ £ which consists 
of the half-circles = {|£| = R, Im ( > 0} and Cf = {\Q = e, lm( > 0} and the intervals 
(s,R) and (-R, —e). We choose R and e such that all of the N negative eigenvalues of H 
lie inside the contour r^ j£ . The function D(( 2 ) is analytic inside and on r^ e as Wj(Q 1S 
analytic in the the upper half-plane Im£ > 0. Thus, 

— DiC 2 ) 

^^ydC = 2niN. (4.14) 

Remember that limim^oo D{C, 2 ) = limi m ^ 00 K(Q/(niC) = 1 + C ( I C I 1 ) - Thus, we can fix 
the branch of the function lnZ)(£ 2 ) by the condition lnl?(£ 2 ) — > as Im£ — > oo. Then, we 
have ln£>(£ 2 ) = ln|L>(C 2 )| + iarg D(Q 2 ). Equation (|4TH|) implies that 

var rKe arg D(( 2 ) = 2ttN. (4.15) 
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We define 77(0) := lim^Q^. rj(k) . This limit exists because of asymptotics (14.61) . It follows 
with (g7rgD that, 

varr H e arg D(( 2 ) = 2(r/(i?) - 77(e)) + var c +arg D(( 2 ) + var c +arg D(( 2 ). 




FIGURE 2. contour of integration Tr^ 



Now, we let R — > oo and e — > 0. Because of (|3.2[) . lim/j^oo var^+ arg L>(C 2 ) = 0. Hence, it 
follows from (|4,15p that 

7/(00) — 77(e) = ttN — -var^+arg D(( 2 ). 
By Proposition 14.51 lim £ _>o var c + arg D(C 2 ) = -(m- 1)tt. Thus, 

( 777- — 1 \ 
iVH ^ — J . (4.16) 

Note that 77(00) = since lnD(( 2 ) — > as \(\ — > 00. It remains to note that in view of 
Theorem 13.21 the following identity is true for A = k 2 , k > 0, 

£(A) = n^ 1 lim arg D(X + is) = ir^ 1 r](X 1 ^ 2 ), A > 0. (4.17) 

£->0+ 

The assertion of Theorem 11.31 then follows by combining (|4.16|) and (|4.17p . 

□ 
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